Thermodynamics of Ising spins on the triangular kagome lattice: Exact analytical 

method and Monte Carlo simulations 
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A new class of two-dimensional magnetic materials CugX2(cpa)e -xEhO (cpa=2-carboxypentonic 
acid; X=F,Cl,Br) was recently fabricated in which Cu sites form a Triangular Kagome Lattice 
(TKL). As the simplest model of geometric frustration in such a system, we study the thermody- 
namics of Ising spins on the TKL using exact analytic methods as well as Monte Carlo simulations. 
We present the free energy, internal energy, specific heat, entropy, sublattice magnetizations, and 
susceptibility. We describe the rich phase diagram of the model as a function of coupling constants, 
temperature, and applied magnetic field. For frustrated interactions in the absence of applied field, 
the ground state is a spin liquid phase with residual entropy per spin so/fcs = | In 72 ~ 0.4752 .... 
In weak applied field, the system maps to the dimer model on a honeycomb lattice, with resid- 
ual entropy 0.0359 per spin and quasi-long-range order with power-law spin-spin correlations that 
should be detectable by neutron scattering. The power-law correlations become exponential at finite 
temperatures, but the correlation length may still be long. 

PACS numbers: 75.30.Kz, 75.40.Mg, 75.10.Hk, 64.60.-i 



I. INTRODUCTION 

Geometrically frustrated spin systems give rise to 
many novel classical and quantum spin liquid phases. 
They may have technological applications as refriger- 
ants, via adiabatic demagnetization^ in which reducing 
the applied magnetic field results in a cooling effect as 
the spins absorb entropy from other degrees of freedom. 
Unlike paramagnetic salts which arc limited by order- 
ing or spin-glass transitions due to residual interactions 
between the spins, geometrically frustrated systems can 
remain in a disordered, cooperative paramagnetic state 
down to the lowest temperatures. Furthermore, they may 
exhibit an enhanced magnetocaloric effect in the vicinity 
of phase transitions at finite applied fields^^^ 

A new class of two-dimensional magnetic materials 
CugX2(cpa)6 • XH2O (cpa=2-carboxypentonic acid, a 
derivative of ascorbic acid; X=F,Cl,Br) was recently 
fabricated, which is an experimental realization of a new 
type of geometrically frustrated lattice. The Cu spins 
in these materials are interconnected in a "triangles-in- 
triangles" kagome pattern, which we refer to as a trian- 
gular kagome lattice (TKL, see Fig. [T]). 

Experiments on the CugX2(cpa)g • XH2O compounds 
show no spontaneous magnetization down to at least 
T = 1.7Kj£ consistent with a spin liquid ground state, in- 
dicating that J aa is antiferromagnetic. However, whether 
J a b is ferromagnetic or antiferromagnetic is still an open 
question. Based on the observation of a partial rather 
than saturated magnetization^ Maruti et al. con- 
cluded that the intcrtrimcr coupling is antiferromagnetic, 
J ab < OJ. In a later theoretical study using variational 
mean field theory on the quantum Heisenberg model on 
the TKL, Strecka£ concluded that the intertrimer cou- 
pling is ferromagnetic, J ao > 0. Regardless, the lack of 
observed hysteresis despite the observation of a magneti- 
zation plateau in finite fielc£ is consistent with a multi- 




FIG. 1: (Color online) The triangular kagome lattice (TKL). 
Solid (open) circles represent "a" ("b") sublattices. Thick 
and thin lines represent interactions J aa and J a b respectively. 
The shaded region represents a unit cell. By shifting the cell 
slightly (as indicated by the dashed parallelogram) it can be 
seen that there are nine spins per unit cell. 



tude of ground states which can be connected by a series 
of local spin flips. 

In this paper we study the classical TKL Ising model 
using exact analytic methods and Monte Carlo simula- 
tions. The purpose of this study is to provide explicit 
predictions, in order to determine to what extent the ex- 
periments can be explained in terms of a classical Ising 
model. In particular, we find significant difference be- 
tween the behavior of the susceptibility as a function of 
temperature for J a b > and J a b < 0, and this may be 
used as an experimental means of distinguishing the two 
cases. Discrepancies between experiments and our theo- 
retical predictions will indicate the vector nature of the 
actual spins (XY or Heisenberg), the effects of quantum 
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fluctuations, or possibly higher order interactions. This 
model was previously studied by one of usi^. In that 
work, Monte Carlo simulations were used to study the 
phase transitions and basic thermodynamics. Zheng and 
Sunii mapped the partition function to that of a kagome 
lattice and found an analytic expression for the phase 
boundary in zero field. 

In this paper we present exact results in zero field at fi- 
nite temperature, and also in an applied field at zero tem- 
perature. We report the full phase diagram as a function 
of coupling constants, temperature, and applied magnetic 
field. We find several field-induced transitions. In partic- 
ular, for frustrated interactions in applied field, we find 
a quasi-long-range ordered phase which maps to hard- 
core dimers on the honeycomb lattice. We complement 
these exact analytic results with Monte Carlo simulations 
on the magnetization and susceptibility We report the 
temperature-dependence of the magnetic susceptibility, 
and show how it can be used to deduce the sign of the 
coupling constants. 

The paper is organized as follows. The TKL Ising 
model is described in Section [TTJ In Section IIII1 wc 
present exact results for the TKL Ising model in zero 
field. In Section [IV] we present exact results at zero tem- 
perature. In Section|V]we present the phase diagram and 
describe the various phases. In Section IVII we present 
Monte Carlo results for the susceptibility, spontaneous 
magnetization, and magnetization curves. In Section IVlII 
we compare our model to models of geometrically frus- 
trated magnets on other lattices, as well as to experi- 
ments on TKL systems. In Section [Villi we present our 
conclusions, and in Appendix we present two mean-field 
approximations in order to illustrate their failure in the 
frustrated regime. 



spins on the a and b sublattices, so that N a : Nb = 2 : 1. 
The space group of the TKL is the same as that of the 
hexagonal lattice, p6m, in Hcrmann-Mauguin notation. 

There are four energy scales in the problem: J aa , J a b, 
T, and h. We have found it most convenient to take \ J a b\ 
as the unit of energy. Thus, the model can be described 
by a three-dimensional phase diagram in the space of the 
three dimcnsionless parameters J a a/\Jab\, h/\J a b\, and 
T/\J a b\- The phase diagram also depends on the sign of 

Jab- 

III. EXACT RESULTS IN ZERO FIELD 

In this section we present exact analytic results for 
the TKL Ising model in zero magnetic field (h = 0), 
including the free energy, internal energy, specific heat, 
and entropy. We use a sequence of V-Y transformations 
and series reductions (shown in Fig. [2]) to transform the 
Ising model on a TKL into one on a honeycomb lattice, 
and then use the known solution for the honeycomb lat- 
tice. (See Ref. [ill for a similar analysis that transforms 
the TKL to a kagome lattice keeping the overall value of 
the partition function unchanged.) We remark that frus- 
trated Ising models with quenched bond disorder may 
be studied by numerical application of V-Y and Y-V 
transformation s 12 i 13 and/or Pfaffian methods^. Our re- 
sults from this section are plotted in Figs. [5] and O using 
the analytic formulas below (as well as series approxima- 
tions for extreme values of T). These thermodynamic 
quantities do not depend on the sign of J a b- 

A. Effective coupling on the equivalent honeycomb 
lattice 



II. MODEL 

The TKL (Fig. [T|) can be obtained by inserting trian- 
gles inside the triangles of the kagome lattice, for which it 
is sometimes referred to in the literature as the "trianglcs- 
in-triangles" kagome lattice. Alternatively, it can be de- 
rived from the triangular lattice by periodically deleting 
seven out of every sixteen lattice sites. This structure 
has two different spin sublattices "a" and "b", which 
correspond to small trimers and large trimers, respec- 
tively. We study Ising spins on the TKL with two kinds 
of nearest-neighbor interactions, the "intratrimer" cou- 
plings J aa and the "intertrimer" couplings J a b- Each 
spin has four nearest neighbors. The Hamiltonian is 



H = -J aa ^2 - Jab a lCT J _ ^ X/ 



(1) 



i£a,j£b 



where <7j = ±1, summations run over the nearest spin 
pairs and all spin sites, h is an external magnetic field. 
The shaded region in Fig. [1] is one unit cell, which con- 
tains 6 a-spins, 3 fe-spins, 6 a-a bonds, and 12 a-b bonds. 
We shall use N a and Nb to denote the total numbers of 



If an Ising model is on a particular lattice contains a 
spin (To connected to only two other spins o\ and oi via 
couplings J\ and J2, then we can 'integrate out' <jq while 
preserving the value of the partition function, in order 
to obtain an effective coupling J12. The transformation 
also produces a constant factor multiplying the partition 
function, A. That is, 



+1 



Ae 



(2) 



for all combinations of values of <j\ and o 2 • It is most 
convenient to write the effective coupling in terms of 
ti,t2,ti2 where U = tanh/3Ji, and the partition function 



changes in terms of x\,X2, £12 where X{ — e 



UU 



-20 Ji 



A = (1 + £1X2) 



X12 
X1X2 



(3) 
(4) 



The x's and t's are related by Mobius duality transfor- 
mations, Xi = (1 — U)/(l + ti) and ti = (1 — Xi)/(1 + Xi). 
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where 



Series 
reductions 



FIG. 2: (Color online). Transformation of a triangular 
kagome lattice (TKL) to a honeycomb lattice. Figure (a) de- 
picts a section of the TKL, with couplings J aa and J a b- The 
procedure begins by applying V-Y transformations to the six 
downward-pointing triangles in each unit cell. This gives lat- 
tice (b). Now take the two strong bonds (J3) in series to give 
J5, and the two weak bonds (J4) in series to give J§, to obtain 
a "3-12 lattice" (c). Apply V-Y transformations to the tri- 
angles to obtain a decorated honeycomb lattice (d). Finally, 
perform series reductions to obtain the honeycomb lattice (f ) 
with a single coupling J;,. 



Similarly, if a spin <tq is connected to only three other 
spins cr 12l 3 via couplings J\2Zi we can 'integrate out' 
Co while preserving the partition function, in order to 
obtain effective couplings J23,3l,i2 together with a free 
energy shift. This is known star-triangle or Y-V 
transformation. The reverse transformation exists, and 
is known as a V-Y transformation: given a "V" of cou- 
plings J23,3i,i2j w e can find an equivalent "Y". Again, 
it is convenient to use the variables t\ = tanh (3J\ (and 
similarly for t 2 , £3) and X\ = e~ 



-2/3 Ji 



y/aia 2 a 3 /a Q 
ti = (cycl.) 



a 1 



(5) 



o 
ai 

A 



1 + ^23^31^12 

£23 + £31^12 (cycl.) 

1 



X1X2X3 



1 + xxxix?, V £232312:12 



(G) 
(7) 

(8) 



and "cycl." means that t 2 , t 3 , a 2 , and 03 are found by 
cyclic permutation of the indices 1,2,3. In general the 
a's may be negative or complex- valued, so that it is not 
correct to replace ^010203/00/01 by y/ 0203/0001. 

Using a sequence of V-Y transformations and series re- 
ductions, we transform the TKL Ising model (with cou- 
plings constants J aa and J a b) into a honeycomb Ising 
model (with a single coupling constant J n ), as shown in 
Fig. O The transformation equations (in terms of the 
ti = tanh (3 Ji variables) are: 



^3 — \J {taa + t a0 2 )/(l + t aa t a b 2 ) 



h 
h 

t 7 

*8 
th 



(tab + t aa tab)/ v/ (t aa + t a t, 2 )(l + t aa t a b 2 ) 



u 2 

l/y/itB+t^ 1 - I) 
t e t 7 

t 8 t 7 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 
(15) 



We can write th directly in terms of t aa and t ao 



ti, 



(1+^aa) tab 2 



(1 - t aa + t aa 2 ) (l + t ab 4 ) - (l - 4t aa + t aa 2 ) t ab 2 



(16) 



It will be convenient to rewrite this in terms of Xi = 



, as this is simpler: 



Xh 



2 (3a; aa 2 + l) x ab 2 



1 



(17) 



For Ising models, both the triangular lattice and the 
kagome lattice can be transformed into the honeycomb 
lattice by a similar procedure. Figure [3] shows the ef- 
fective dimensionless coupling of the honeycomb lattice, 
(3Jh, as a function of the inverse temperature for antifer- 
romagnctic Ising models on different lattices for compar- 
ison. Note that: 

• For the triangular Ising AF (TIAF), the effective 
honeycomb coupling is imaginary, and there is no 
long-range order. 

• For the kagome Ising AF (KIAF), the effective hon- 
eycomb coupling is real. As the kagome couplings 
are increased, the effective honeycomb coupling in- 
creases until it saturates at large /3. However, it 
never grows beyond the critical coupling for the 
honeycomb model, /3Jf = tanh _1 l/\/3 w 0.658. 
Therefore the KIAF does not have a phase transi- 
tion. 
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FIG. 3: (Color online). Effective dimensionless coupling f3Jh 
of the honeycomb lattice Ising model as a function of the 
original inverse temperature j3, for the triangular Ising anti- 
ferromagnet (TIAF), kagomc Ising antiferromagnet (KIAF), 
or TKL Ising antiferromagnet (TKLIAF). In the case of the 
TIAF, (5Jh is imaginary, and Im (5Jh is shown as a dashed 
curve; for the KIAF and TKLIAF, /3J h is real. For the TK- 
LIAF in the unfrustrated regime, j3J n crosses the dotted line 
(the critical coupling of the honeycomb Ising model) , indicat- 
ing a phase transition. The other models do not have phase 
transitions. 
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FIG. 4: (Color online) Phase diagram of the TKL Ising model 
in the (J a a,T) plane, for J a t = ±1 and h = 0. The thick 
curve is the exact solution (equivalent to that in Ref. ITlh . 
The dotted and dashed curves are mean-field approximations 
(see Appendix, Sec. [A}. The ordered phase is ferromagnetic 
if Jab > and ferrimagnetic if J a b < 0. The disordered state 
is paramagnetic. For J aa < —J a b, it persists down to T — 0, 
where the entropy is In 72 per unit cell and the susceptibility 
obeys a Curie law. The arrows indicate the two TKLIAF 
cases discussed in Fig. [3] 



1/Pc of the TKL Ising model: 



• For the triangular kagome lattice Ising antiferro- 
magnet (TKLIAF) in the unfrustrated regime (e.g., 
Jaa = — 0.8, J ao = —1), the plot of (JJh intersects 
the dotted line, indicating a phase transition at 
/?«3.5. 

• In contrast, for the TKLIAF in the frustrated 
regime (J aa = -0.8, J & = -1), f3J h grows to a 
maximum and then decays to zero, indicating the 
absence of a phase transition. Paradoxically, and 
unlike the other lattices, stronger bare couplings 
lead to a weaker effective coupling! The fact that 
(iJh — > as (3 — > oo is closely connected to the fact 
that the residual entropy of the TKL lattice has 
the simple value of In 72 per unit cell (as we will 
show), unlike the cases of the TIAF±£ and KIAF±£ 
in which the residual entropies are non-trivial two- 
dimensional integrals. 

B. Phase boundary 

The phase boundary of the TKL Ising model in zero 
applied field was calculated in Ref. [111. We show an alter- 
native exact derivation of the results here for pedagogical 
reasons. Once we have used the techniques of Sec. IIII Al 
to map the TKL Ising model into the Ising model on a 
honeycomb lattice, we can use known results for the hon- 
eycomb lattice. The critical temperature of the honey- 
comb Ising model is given by t c h = tanh/3J^ = l/v3, or, 
cquivalently, x c h = cxp—2[3J h = 2 — V3. Substituting in 
the equivalent coupling of the honeycomb lattice, Eq. ll7[ 
leads to an implicit equation for the critical temperature 



e -4/3 c J ao = (a /3 - 1) cosh4/? c J ab - (VS + l) . (18) 

Eqn. 7 of Ref. [ll| is equivalent to the simpler ex- 
pression reported here. This critical curve is plot- 
ted in Fig. [H For large ferromagnetic J aa , the criti- 
cal temperature saturates at a finite value, T c /\J ao \ sa 

4/ In (2 + V3 + VG + 4V3") w 2.00838. As J aa is re- 
duced towards — 1, the critical temperature falls to zero. 
Near J aa = — \J a b\, the critical curve is approximately 
linear: 

Since the mapping from the TKL to the honeycomb 
lattice preserves the nature of the singularity in the parti- 
tion function, the phase transition is a continuous second- 
order transition in the 2D Ising universality class. 

In zero field (h = 0) the partition function Z is invari- 
ant under a change of sign of J a0l and the topology of 
the phase diagram is independent of this sign, although 
the identification of the phases is not. First consider 
the model for the case J aa = 0, where the TKL reduces 
to a decorated kagome lattice. If J a b is ferromagnetic, 
the model develops ferromagnetic order below the Curie 
temperature. If J ao is antiferromagnetic, the model de- 
velops ferrimagnetic order below the ordering tempera- 
ture. Although the decorated kagome lattice is bipartite, 
the numbers of spins on the a- and 6-sublattices are not 
equal: there are twice as many a-spins as 6-spins. Hence 
we have ferrimagnetic order, with unequal numbers of 
up and down spins producing a net moment. However, 
the transition temperature, free ener:gy, internal energy, 
specific heat, and entropy are independent of the sign of 
J a b in the absence of applied field. Now introduce the 
coupling J aa . If this is ferromagnetic, it has very little 
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effect, since the a-spins already have a tendency to align. 
However, if J aa is antiferromagnetic, it fights against the 
ordering induced by J a b- If Jaa is sufficiently antiferro- 
magnetic, (J a a/\Jab\ < ~ 1); the system is in a frustrated 
regime with no order even at zero temperature. 



In order to get accurate numerical results one has to fur- 
ther split the range of integration according to the sin- 
gularities of the integrand. 

Thus the partition function of the TKL Ising model 
(per unit cell) is 



C. Partition function 

The partition function per unit cell, ztkl> of the TKL 
Ising model is equal to that of the equivalent honey- 
comb lattice zh, multiplied by the factors Zi, . . . ,zq be- 
low which are accumulated during the sequence of \7-Y 
transformations and series reductions necessary to effect 
the transformation. These factors are 



(first V - Y) 




1 + X$X7 V Xg,X7 



(19) 
(20) 
(21) 

(22) 
(23) 
(24) 



The total accumulated partition function change is there- 
fore 



ZTKL = Zi 6 Z 2 6 Z 3 3 Z 4 2 Z 5 3 Z 6 3 Z H 



^3 



(1 + x ab 2 ) 2 (1 + 6x aa 2 x ab 2 + x a b A ) 

l + 2{l + 6Xaa 2 )x ab 2 +X ab 4 



ZH 



(25) 



The partition function per unit cell of the honeycomb 
lattice has been calculated in the literature by, e.g., the 
Pfaffian method^. It is 

z H (x h ) = ^ l - ' Xh ^ exp [%Sl{w{x h ))] (26) 

where 

, f 2 " dp f 27r dq. , 
\l[wj = / — / — m(w — cosp — cos q — cos(p + qn 
J a 27r J 2tt 

(27) 

and 



w(x h ) 



1 - 2x h 



Qx h 2 - 2x h 3 



x h 



2xk{^ - x h ) 2 

For the purposes of numerical evaluation, we rewrite 
the function Q(w) in the following form 



72 



dpln 



cosp + arccosh 



w — cos 2p 
2 cosp 



(29) 



ZTKL 



{Xaa,X a b) = * (x a a , X ab ) exp [^tt(w(x h {x aa , X a b) ) )] 

(30) 



where 



* = 2x aa 6 x ab (1 - Xab 1 ) x 

+ 3x aa 2 ) (1 + Gx aa 2 x ab 2 + x ab A ) 



(31) 



andfi, w, andx h are defined in Eqs. f(|27]l). ((EH)), (IffUl), 
respectively The total partition function Ztkl is related 
to the partition function per unit cell ztkl by .Ztkl = 
z tkl ' where N is the number of unit cells. We show plots 
of —InzTKL /{f3\J aa \) in Figs. [5] and [5] (red curves) in the 
unfrustrated and frustrated regimes, respectively. 



D. Energy 

The energy per unit cell of the TKL Ising model can 
be obtained by differentiation of the partition function: 



dlnz 
dfT = ' 

^ J %Xi 

i—aa^ab 



dx aa d In z dx a b d In z 

dP dx aa d/3 dXab 

d\n*i> dxh dw dil 



dxi 



dxi dxh dw 



(32) 
(33) 



where \1/ is given in Eqn. [31] 1^ is the Green function of 
a particle on a triangular lattice and can be expressed in 
terms of the complete elliptic integral of the first kind, 



dn_ f 27T dp r** dq 

dw J 2tt J 2tt w — cosp — cos q — cos(p + q) 



1 



"tt(-io-1) 3 / 4 (-m>+3) 1 / 4 ' 



(34) 
(35) 



L \2 + 2(w+l)(-w-l) 1 ^ 2 (-w+3) 1 / 2 J • ( 36 ) 



(28) We show plots of —u/|J QQ | in Figs. [5]and[n](green curves). 



E. Specific heat 

The heat capacity per unit cell, c 
taincd by further differentiation: 



can be ob- 



6 



2{ (PJaa) 2 X a 



d 2 ln^ ( 8x h d 2 x h 



9x aa 0x aa \dx a 



dx n 



'O' 



dx h 

dXaa 



x aa (w"Q' + (w') 2 n") 



+ [previous term with aa replaced by ab] 
d 2 In # d 2 x h 



2/3 JaaJab%aa%ab 



dx aa dx a i, 9x aa 0x a i) 



dx h dx h 



( w "n' + {w') 2 n") 



(37) 



where w' , CI', etc., represent derivatives of the functions 
w{xh) and il(w) with respect to their arguments. Wc 
show plots of c in Figs. [5] and [6] (blue curves). 

In the unfrustrated case (J aa > 0, Fig.O, the specific 
heat has a broad hump just above T — J aa , and a sharp 
peak near T = 2\J a b\- The broad hump is due to ferro- 
magnetic alignment within each a-plaquette. The sharp 
peak corresponds to the phase transition to a ferromag- 
netic (for J a b > 0) or ferrimagnetic (for J ao < 0) state, 
governed by the weakest links, J a b- The position of the 
sharp peak is consistent with the solution of Eqn. [18] 
(See also Fig. [4]) In the frustrated case, broadened fea- 
tures remain at both of these energy scales, as shown in 
Fig. E 



F. Zero-temperature limit: residual entropy 

Far from the critical curve, results for two limits can 
be obtained, corresponding to the ordered phase and to 
the disordered phase (which persists even at zero tem- 
perature). 

In the first case, J a a/\Jab\ > — 1 ; the system orders 
at low temperatures, going into cither the ferromagnetic 
state (for J a b > 0) or ferrimagnetic state (for J a b < 0). 
In the low temperature limit, the partition function and 
internal energy may be expanded as 



In Z{j3) = 6 In x a b - 3 In x aa + 



6x 2 



(38) 



u{p) = -12\J ab \ + 6\J aa \ + 24e^ J ^- J ^\ J ab - J aa \ - 

(39) 

As T — > 0, the residual entropy is zero, whether in the 
ferromagnetic phase or the ferrimagnetic phase. 

Suppose J aa /\Jab\ < — 1- In this case, the model 
becomes frustrated when T — > 0, (i — » oo. We have 
x aa,Xab °o> x h — * 1 — 5 w — > oo. Expanding Cl(w) w 
In w in a series in w, and then expanding In Z as a series 
in x aa and x a b, we find: 



In Z(/3) = In 72 + In x aa + 2 In x a b 



1 



%ab 



%ab 



3x aa &%aa' 



+ 



(40) 



proximation for the energy per unit cell: 



u(0) = -2\J aa \ - 4|J ah | + | e -^U«-J-l| J aa - J ab \ 



-40\Jaa +Jab\ 



Jab 



i„-4/3|J a 
3 e 



I Jaa | 



+ 8e -4/3|J Qb || Jah | + 



(41) 



The first two terms are the ground state energy. The 
other terms represent different types of excitatiotns 
about the ground state (s). The coefficients represent rel- 
ative degeneracies of excited states and ground states, 
and the exponents represent excitation energies. 

The first term in the series expansion of In Z gives the 
residual entropy per unit cell: 



so 



lim lim QnZ+pu) = In 72. (42) 



Thus the residual entropy is exactly In 72 ss 4.2767 . . . 
per unit cell, or ^ In 72 w 0.4752 . . . per site. This num- 
ber will be discussed in more detail in Sec. IIVI 




FIG. 5: (Color online). Thermodynamic functions vs temper- 
ature T for unfrustrated (ferromagnetic) couplings J aa — 5 
and | J a j,| = 1. Red: free energy f(T) = — fcsTlnZ. Green: 
energy u(T). Blue: heat capacity c(T). Yellow: entropy s(T). 
All values quoted per unit cell; each unit cell contains 9 sites. 



We can thus obtain the following low-temperature ap- 
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^ T 



5 10 15 20 



FIG. 6: (Color online). Thermodynamic functions vs temper- 
ature T in frustrated regime, J aa = —5 and \J a b\ = 1, where 
the intratrimer coupling J aa is strong and antiferromagnetic. 



IV. EXACT RESULTS AT ZERO 
TEMPERATURE 



urations per unit cell to the ground state manifold. For 
any given configuration of the 6-spins, each a-trimer in 
the lattice has a 3-fold degeneracy. Since there are two 
a-trimers per unit cell, these contribute a factor of 3 2 = 9 
to the ground state degeneracy. The total degeneracy per 
unit cell in the ground state is therefore 8 x 9 = 72, as 
we showed in Sec. IIIIFI 

The fact that the 6-spins arc effectively independent 
also means that the correlation function is "perfectly lo- 
calized" : it is exactly zero beyond a distance r b (the 
distance between two 6-spins). The correlation length £ 
is thus zero (where £ is defined as the asymptotic decay 
length of the correlation function at large distances). 

For comparison, at T = h = 0, the triangular Ising 
AF has power-law correlations. The kagome Ising AF is 
more frustrated than the triangular lattice case, since its 
ground state has exponentially decaying correlations. We 
have shown here that the ground state of the TKL Ising 
AF in the frustrated regime has perfectly localized cor- 
relations, making this model even more frustrated than 
cither the triangular or kagome cases. 



In this section we show how the zero-temperature 
phase diagram (along with the thermodynamic proper- 
ties and correlations of the various phases) can be sys- 
tematically deduced with and without applied field by 
considering ground states of large triangular plaquettcs. 
By explicitly comparing the ground state energies, we de- 
rive the phase diagram. The results of this section are 
summarized in Figs. IT01 and [Til 



A. Zero Field (phases V and VI) 

The phase diagram for zero applied field is shown in 
Fig. [U For J aa > — \J a b\ and low temperature, the sys- 
tem is in an ordered phase which is ferromagnetic for 
J a b > 0, and ferrimagnetic if J a b < 0. For J aa < — \J a b\, 
the system remains disordered even at zero temperature, 
with a residual entropy of sq = ln72. The zero-field dis- 
ordered phase is labeled phase V in Figs. [TD1 and [TTl and 
the zero-field ordered phase is labeled phase VI. 

The degeneracy of the ground state manifold can be 
understood by considering the energetics of a single large 
plaquette, i.e., an a-spin trimer along with its enclosing 
6-spin trimer. Representative plaquette configurations 
within the ground state are shown in Fig. [7] We enu- 
merate all possible plaquette energies in Table |TJ As can 
be seen from the table, for any of the 2 3 = 8 possible 
configurations of the 6-trimer, there are three and only 
three configurations of the enclosed a-trimer which are all 
within the ground stated This means that the 6-spins 
are effectively free within the ground state manifold. 

In order to count the ground state degeneracy, we now 
turn to the unit cell. Since the 6-spins are free, and there 
are 3 6-spins per unit cell, this contributes 2 3 = 8 config- 
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TABLE I: Energy of a large triangle (consisting of 3 6-spins 
and 3 a-spins) as a function of the configurations of a and b 
spins. For clarity of presentation, we present the results for 
the case J aa = —3, J a b — — 1, h = 0, but the form of the 
table is representative of the entire line J aa /\ J a b\ < — L The 
boldfaced numbers indicate the lowest-energy configurations, 
which have energy J aa + U a b- Note that each row contains 
exactly three boldfaced numbers. See text for discussion. 



B. Saturated ferromagnetic phase (phase I) 

At very high fields, when h/\J a b\ > 
max(4, 2| J aa |/| J a b| +2), we find that there is a unique 
ground state where all the 6-spins and a-spins are up. 
This state is easily seen to have magnetization m = 9, 
entropy s = 0, and energy u = —3h + 6 |J aa | + 12 \ J a b\ 
per unit cell. 
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(a)All three ft-spins pointing up 
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(b)Two of three b-spins pointing up 




FIG. 7: (Color online). The ln72 phase. Given a configu- 
ration of the three 6-spins (on the outer triangle), it can be 
shown, by enumerating all possibilities, that there are exactly 
three states of the a-spins (on the inner triangle) . The figure 
illustrates this for two configurations of 6-spins; results for the 
other configurations can be seen from the Ising symmetry and 
the local rotational symmetry of the triangular plaquettes. 



C. Ferrimagnetic phase (phase II) 

At lower fields < h/\J a b\ < 4 and when J aa /\J a b\ > 
— 1, there is a ferrimagnetic phase: the a-spins, being 
more numerous, align parallel to the field (up); the b- 
spins are then induced to point down due to the antifer- 
romagnctic J a \, interaction. This phase has m = 3, s = 0, 
and u = —2>h + 6 | J aa \ — 12 \ J a b\ per unit cell. 



D. Log 9 phase (phase III) 

If 4 < h/\ J a b\ < 2| J aa \/\ Jab\ + 2, the 6-spins are com- 
pletely polarized, but each a-triangle has 3 degenerate 
states (see Fig. [5J . Therefore, in this phase the system 
is equivalent to a set of non-interacting three-state Potts 
spins. This phase has m = 5, s = ln9 = 2.1972 . . . , and 
u = —5h — 2 | J aa \ +4 | J a f, | per unit cell. Again, the corre- 
lation function is perfectly localized, and the correlation 
length is £ = 0. 



E. Dimer phase (phase IV) 
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FIG. 8: (Color online). Phase III of the TKL Ising model. 
When the field is quite strong, the b-spins (outer triangle) 
are all polarized, and each a-trimer (inner triangle) has three 
degenerate ground states. 



like a lattice gas on the kagome lattice with nearest- 
neighbor exclusion, at maximal density. The ground 
states can be mapped to configurations of dimers occu- 
pying a honeycomb lattice (see Fig. [5]). This problem 
has been solved exactly using the Pfaffian method^; the 
entropy per unit cell is 



— - / dp dq In (l — 4 cosp cos q + 4 cos 2 g) 
8tt 2 Jo Jo 



0.3231 



(43) 



Therefore, the entropy of the TKL phase IV is also 
0.3231 . . . per unit cell, or 0.03590 . . . per site. 



A A & 

TKL phase IV 




Honeycomb dimer model 



FIG. 9: (Color online). Phase IV of the TKL Ising model: 
Mapping of the spin configurations to configurations of dimers 
on a honeycomb lattice. The yellow patches represent defects 
in the spin configuration, which correspond to vacancies (or 
monomers) in the dimer picture. 



The most interesting situation occurs when < 
h/\J a b\ < 4 and J a a/\Jab\ < — 1. Table U shows that the 
system will have the lowest energy if each &-trimer has 
exactly one &-spin pointing down. Counting the number 
of ways to satisfy this constraint globally is a non-trivial 
problem. The situation is the same as that for a kagome 
Ising AF at T = and < h < 4|Jj ia g ome | ) which has 
been studied befor o 20 ' 21 i 22 ' 23 . The down-6-spins behave 



The correlation function of the 6-spins, Cfcb(r) = 
(er{,(0)<7(,(r)}, is equivalent to the dimer-dimer correlation 
function. This correlation function decays as a power law, 
1/r 2 , and nice visualizations in real space and reciprocal 
space are given in Ref. Thus, the model is in a criti- 
cal phase, in the "Villain-Stcphenson universality class" . 
The magnetization per unit cell is to = 3 and the energy 
per unit cell is u = —3h — 2 \J aa \ — 4 \J a b\- 
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FIG. 10: (Color online). Phase diagram of the TKL Ising 
model in the (J aa ,h) plane, for J a t, < (antiferromagnetic 
intertrimer coupling) and T — 0. The phase diagram is sym- 
metric under a change of sign of h. The phases P, IP, IIP, IV 
are just mirror images of I, II, III, IV obtained by swapping 
up and down spins. The yellow region represents the critical 
phases with power-law correlations; note that it does not in- 
clude the thin line at h = 0, which is phase V, the In 72 phase 
described in the text. The thick line (VI) persists as a true 
ferrimagnetic phase transition at finite T; all other lines turn 
into crossovers. 



F. Phase diagram 



We now combine the above results in order to re- 
port the full phase diagram of the TKL Ising model. 
Fig. [10] shows the phase diagram for antiferromagnetic 
intertrimer coupling, J a b < 0, and in Fig.[TTJwe show the 
phase diagram for ferromagnetic intertrimer coupling, 
J a b > 0. In both cases, the phase diagram is symmet- 
ric under for h — ► —h, with simultaneous change of sign 
of all spins. The entropy and magnetization change dis- 
continuously across every zero-temperature phase bound- 
ary. Exactly on the phase boundaries and intersections 
of these boundaries, the entropy will be higher than in ei- 
ther adjacent phase, because the system can choose from 
states within each set of ground states. 

Note that the high-field phases (I and III) are common 
to both phase diagrams. The more interesting case is 
that of Fig. I10[ which has antiferromagnetic intertrimer 
coupling J ao < 0, and more phases at intermediate field 
strength. Right at h = in both phase diagrams, the 
ground state is phase V, the spin liquid with residual 
entropy so = ln72 that we discussed in Sec. IIV Al When 
Jab < 0, the application of an infinitesimal field induces a 
critical state with power law correlation functions (phase 
IV), which we have mapped to the problem of hard-core 
dimers on a honeycomb lattice. 
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FIG. 11: (Color online). Phase diagram of the TKL Ising 
model in the (J aa ,h) plane, for J a t > (ferromagnetic inter- 
trimer coupling) and T = 0. 



V. FINITE TEMPERATURE AND FINITE 
FIELD 

Most of the phase transitions in the zero-temperature 
phase diagrams are destroyed by thermal fluctuations. 
The clear exception is phase VI, which has spontaneously 
broken Z2 symmetry. This long-range ordered phase sur- 
vives at finite temperature, and has a true phase transi- 
tion at a Curie temperature T c ^ 0. Since h is a relevant 
perturbation, this finite-temperature phase transition is 
destroyed at any finite field, leaving only a crossover. 

The other transitions in Figs. [TU1 and [HI arc not char- 
acterized by a spontaneously broken symmetry with an 
order parameter. Therefore, they cannot persist at fi- 
nite T as traditional order-disorder transitions. How- 
ever, a more subtle analysis is required to understand 
whether the critical phase, phase IV, persists at finite 
T (bounded, e.g., by a Kosterlitz-Thouless transition 
curve). The situation is quite similar to that for the 
kagome Ising AF described in Ref. [23|. Based on the 
table of ground state energies, Table [TJ it is possible to 
enumerate the types of defects that can occur in phase 
IV. These defects correspond to breaking a dimer into 
two monomers (see Fig. [S]), and can only be created in 
pairs. Each defect has energy AU = mm(h,A\J a b\ — h). 
The entropy AS associated with creating a defect pair is 
at least 0(lnL 2 ), because one may choose to break any 
of the 0(L 2 ) dimers; in fact the entropy is even greater 
than this because the resulting monomers can be moved 
apart, resulting in many new configurations. Hence at 
any finite T, the density of monomers is finite. A the- 
orem of Licb and Hcilman 24 states that monomer-dimcr 
models cannot have phase transitions at any finite density 
of monomers. Therefore, the critical phase only exists at 
T = 0, and it is destroyed at finite T . 

However, if the density of monomers is low, the corre- 
lation length £ may still be long; naively, one might ex- 
pect £ 2 - l/n monomcrs cx exp(£ , monomor /T), but a more 
careful treatment accounting for the effective Coulomb 
attraction between monomers gives the result £ 2 ~ 
l/«monomcrs oc exp(8£ monomcr /7T) & This scenario is 
consistent with the size-dependent peaks in the susccp- 
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FIG. 12: (Color online) .Temperature dependence of heat ca- 
pacity and susceptibility from Monte Carlo simulation for 
Jaa = 5 and J a b = +1 with system sizes L = 12,24,36,60. 
Red line (no dots) represents the exact solution of specific 
heat. Peaks become taller and narrower as L increases. 



tibility that we find from Monte Carlo simulations (see 
Sec. ED). 
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VI. MONTE CARLO SIMULATIONS 



In this section we present the results of Monte Carlo 
(MC) simulations of the TKL Ising model, Eq. (fTJ) , for 
various combinations of parameters. The simulations 
corroborate our analytic predictions and also allow us 
to perform calculations at finite h and T as well as to 
compute the magnetization and susceptibility. 

We use the Wolff algorithm for J aa > at h = 0, and 
the Metropolis algorithm for J aa < at h = and h ^ 0. 
The system sizes we use are L = 12, 24, 36, 60, where L is 
the length of the underlying triangular lattice, so that the 
total number of spins is N = j^L 2 ; for periodic boundary 
conditions, L should be a multiple of 3 in order to avoid 
introducing boundary defects and additional frustration. 

In order to evaluate the heat capacity C and magnetic 
susceptibility \, we use the fluctuation-dissipation theo- 



FIG. 13: (Color online). Temperature dependence of heat 
capacity and inverse susceptibility (l/x(T)) and T\(T) from 
Monte Carlo simulation for J aa = — 10 and J a b = — 1 with 
system sizes L — 12, 24, 36, 60. 
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where (H) and (M) are the Monte Carlo averages of the 
total energy (i.e., the Hamiltonian) and magnetization, 
respectively. We define the sublattice magnetizations as 



a - N a 



(46) 
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FIG. 14: (Color online). Temperature dependence of heat ca- 
pacity, inverse susceptibility (1/%(T)) and T\{T) from Monte 
Carlo simulation for J aa = — 10 and J a i — 1 with system sizes 
L = 12,24,36,60. 

mb = au ' ( 47 ) 

where 7V a is the number of a-spin sites, and ./V;, is the 
number of 6-spin sites. 

A. Zero magnetic field 

We first show Monte Carlo results at zero field. Fig.fTSl 
shows the temperature evolution of the heat capacity and 
susceptibility in a representative unfrustratcd case, J aa — 
5J a b > 0. The Monte Carlo results for the heat capacity 



are consistent with the exact results in Sec. IIVI the peak 
in the heat capacity in our simulations becomes taller 
and narrower as L increases, tending towards the exact 
solution for L = oo. 

Figs. [12] and [JJ] show results for two frustrated pa- 
rameter combinations (J a a/\Jab\ < — !)• In both cases, 
the susceptibility shows a marked difference from the 
case of ferromagnetic intratrimcr coupling, J aa > 0. For 
Jaa > 0, the susceptibility shows a sharp peak at the 
Curie temperature T c , whereas for J aa < 0, it tends to 
oo as T -> 0. 

Unlike the free energy, heat capacity, entropy, and in- 
ternal energy, the susceptibility depends on the sign of 
J ao . Our predictions about the susceptibility can be used 
to distinguish whether a physical TKL Ising system has 
ferromagnetic or antiferromagnetic coupling between the 
a and b sublattices. For J a b > 0, the inverse susceptibility 
l/x(T) shows two linear pieces, with a crossover at the 
lower coupling constant J a b- However, when J ao < 0, the 
inverse susceptibility approaches zero much slower than 
the J a b > case. 

The difference between antiferromagnetic and ferro- 
magnetic coupling J a b is even more evident in the plots 
of T\- At low temperatures, T\ tends to a small finite 
constant as T — ► if the intertrimer coupling is antifer- 
romagnetic, J a b < 0; whereas, if the intertrimer coupling 
is ferromagnetic, J a b > 0, Tx goes through a minimum 
around T ~ J Qa , increases sharply, and saturates at a 
finite value as T — > 0. 



B. Finite magnetic field 

In Figs. [T5] and [TBI we show Monte Carlo results at 
finite field h and low temperature T = 0.1 for J aa = —2 
and J a b = — 1- The magnetization curves in Fig. 1151 have 
a series of steps and plateaux (as is typical of frustrated 
spin systems). Starting from h = 0, as a field is applied, 
both the a and b sublattices immediately respond as the 
critical dimer phase is induced, developing (normalized) 
sublatticc magnetizations of m a = m = 1/3. As field 
is increased, the 6-spins are more easily polarized, while 
the a-sublattice spins only become fully polarized when 
the magnetic field is strong enough. 

The susceptibility as a function of applied field has 
a series of peaks, as shown in Fig. [TH] The peaks at 
h = ±4 increase with increasing system size from L = 12 
to L = 24. This indicates that the correlation length 
(in the vicinity of the peaks) is comparable to, or larger 
than, the system size. This does not indicate a true finite- 
temperature phase transition. Rather, because phase IV 
is critical (with infinite correlation length), the correla- 
tion length diverges as temperature is lowered toward 
this phase. At the low temperatures we have simulated, 
the correlation length is comparable to our system sizes. 
See Sec. El 
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FIG. 15: (Color online). Magnetization m vs. external field 
from Monte Carlo simulation for J aa = — 2, J a j, = —1 at 
T — 0.1 for L — 12. Sublattice magnetizations m a and 
are also shown. 
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FIG. 16: (Color online). Susceptibility x vs - external field 
h from Monte Carlo simulation for J aa = —2, J a b = — 1 at 
T = 0.1 for L = 12. 



VII. DISCUSSION 

Phases I, II, and III can easily be simulated using the 
Metropolis algorithm. However, phase IV has extremely 
slow dynamics when simulated using the Metropolis al- 
gorithm at low temperature - this could be described 
as glassy dynamics. We believe that the simulations 
may be made more efficient by using geometric cluster 
methods 2 ^ or by augmenting the Metropolis algorithm 
with directed-loop updates tailored to the honeycomb 
dimer state i 26 ! 27 The slow dynamics under the Metropo- 
lis algorithm may, however, be representative of the true 
dynamics in a physical realization of the TKL Ising an- 
tiferromagnet. 

It may be possible to stabilize more phases at finite 
temperature by introducing an appropriate perturbation. 
For example, introducing coupling in the third direction 
may be sufficient to stabilize the critical phase it finite 
T. In addition, it may be possible to induce a Kasteleyn 



transition in phase IV by applying some kind of orienting 
field {e.g., uniaxial strain) to the TKL, similar to what 
was suggested by Moessner and Sondhi 2 ^ on the Kagomc 
lattice. In light of existing studies on triangular and 
kagome lattices, it is reasonable to expect that adding 
next-nearest-neighbor interactions to the TKL will pro- 
duce an even richer phase diagram, like Kagome latticed 



Comparison with Other Frustrated Ising 
Models 



In Table [TTJ we present a comparison of the spin-spin 
correlations and the residual entropy among the frus- 
trated TKL, triangular, and kagome lattices. The sim- 
plest example of a geometrically frustrated system is the 
triangular Ising antiferromagnct (TIAF), i.e., a set of 
Ising spins on a triangular lattice with antifcrromagnctic 
pairwisc couplings. Due to the presence of odd cycles in 
the lattice graph, it is impossible for all pairs of nearest 
neighbor spins to be simultaneously antiparallcl. As a 
result, the antiferromagnetic interactions are unable to 
produce long-range order even at zero temperature. In- 
stead, at zero temperature, the TIAF has a quasi-long- 
rangc-ordercd state in which the correlations decay with 
distance as j-- 1 / 2 29 ^ This ground state is macroscop- 
ically degenerate, with a zero-temperature residual en- 
tropy of 0.3231fc_B per spin^i^L^^ This number is the 
same as the entropy per unit cell of the random dimer 
model on a honeycomb lattice, and it crops up in many 
other places. Although the TIAF has no phase transition 
in zero field, the application of a finite field produces a 
surprisingly rich phase diagram. A small field induces 
a Bcrczinskii-Kostcrlitz-Thoulcss (BKT) transition to a 
'spin crystal' phase that breaks translational symmetry 
and has long-range order - the correlation function oscil- 
lates with distance but does not decay. At larger fields 
the crystalline order is destroyed via a transition in the 
3-state Potts universality class. In this limit the TIAF is 
related to Baxter's exactly soluble hard hexagon lattice 
gas models 

Another frustrated spin system is the antiferromag- 
netic Ising model on a kagome lattice, formed by peri- 
odic removal of a quarter of the sites from the triangular 
lattice. Unlike the TIAF, the KIAF in the absence of 
field has pair correlations that decay exponentially at all 
temperatures including T = 0i 35 i 36 Its ground state en- 
tropy is 1.5055&B per unit cell or 0.5018/cb per spin^ 
higher than that of the TIAF because the smaller coor- 
dination number allows more freedom. At finite h, there 
is a different spin liquid state that can be mapped to 
random dimers on a honeycomb lattice. In this state 
the spin-spin correlation function decays as a power law, 
1/r 2 . The residual entropy is 0.3231fcs per unit cell, or 
0.1077fc B per spin. 

The ground state of the TKL Ising AF without applied 
field is even more frustrated than that of the kagome lat- 
tice: the correlation function becomes exactly zero be- 
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long-range-ordered 
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In 72 


for r > r bb 


TKL in weak field 
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TABLE II: Comparison of various frustrated Ising models at 
T = 0. Residual entropies are quoted per unit cell; they are 
the logarithms of irrational numbers, unless otherwise stated. 



yond a certain cutoff radius. The residual entropy per 
spin is ^ In 72 = 0.4752 .... At finite h, there is a corre- 
lated, critical spin liquid state which we have mapped to 
hard core dimers on a honeycomb lattice. This state has 
a residual entropy of 0.3231fcs per unit cell, or 0.03590 
per spin. (See Table llll) 



B. Comparison with Experiment 

As discussed in the introduction, the recently fabri- 
cated family of compounds CugX2(cpa)g • XH2O have an 
arrangement of the copper sites which forms a triangular 
kagome lattice. There is no evidence of spontaneous mag- 
netization down to at least T = 1.7K*£ consistent with a 
spin liquid ground state, which indicates that J aa is an- 
tifcrromagnctic. However, there is no agreement yet as 
to whether J a b is ferromagnetic or antiferromagneticj 7 -*^ 

In studying the magnetic susceptibility, several groups 
find that the slope of 1/x versus T is roughly linear 
at high temperature, but as T is lowered, the slope 
increases . 7 ' 8 ! 37 . This is consistent with cither J a b fer- 
romagnctic or antiferromagnetic, as seen in Figs. |13(b) 
and |14(b)| A distinction can be made, though, by the be- 
havior of Tx- Whereas T\ saturates to a finite value at 
low T for antiferromagnetic intertrimer coupling J a b < 0, 
Tx reaches a minimum at intermediate T before saturat- 
ing at a finite value as T —> if the intertrimer cou- 
pling is ferromagnetic, J a b > 0. To the extent that the 
CugX2(cpa)6 • XH2O materials can be described by an 
Ising TKL model like the one in this paper, our calcula- 
tions indicate that the intertrimer coupling must be fer- 
romagnetic, J a b > 0, so that upon application of a field, 
the system should be in Phase III, rather than Phase 
IV. Since there is a striking difference in the residual en- 
tropies in these two phases (sq = In 9 = 2.1972- •■ in 
Phase III vs. so = 0.3231 • • • in Phase IV), heat ca- 
pacity measurements will also be useful in distinguishing 
these phases. Other future experiments, including neu- 
tron scattering, NMR, and /iSR, can also provide useful 
data for comparing to models of frustrated magnetism on 
the TKL. 



VIII. CONCLUSIONS 

In conclusion, we have studied an Ising model on 
the triangular kagome lattice (TKL) with two interac- 
tions J aa and J a b, temperature T, and external field 
h. We have calculated the complete phase diagram in 
{J aa , Jab, h, T) parameter space exactly. Furthermore, we 
have obtained exact results for thermodynamic quantities 
(free energy, energy, heat capacity, and entropy) at all T 
for h — 0, and at all h for T = 0, and plotted them for 
representative cases. 

In the experimentally relevant regime, |J aa | 3> | «7a.6 1 , 
if J aa is ferromagnetic, the specific heat shows a broad 
hump corresponding to intra-trimer ordering, as well as 
a sharp peak at lower temperatures due to the onset of 
true long-range order. If J aa is antiferromagnetic there 
are two broad features. 

We have computed the magnetization M(T, h) and sus- 
ceptibility x(T, h) in various regimes using Monte Carlo 
simulations. To the extent that experiments on the 
CugX2(cpa)g • XH2O materials can be compared with an 
Ising TKL model, our calculations indicate that J ao > 0. 

We find that at strong frustration and zero field, as 
temperature is reduced, the model enters a spin liquid 
phase with residual entropy sq = In 72 per unit cell, with 
"perfectly localized" correlations. This stands in contrast 
with the triangular and kagome Ising antiferromagnets, 
whose residual entropies cannot be expressed in closed 
form. 

The most interesting feature of the model is a corre- 
lated critical spin liquid phase (with power-law correla- 
tions) that appears at strong frustration, weak fields, and 
zero temperature. We have mapped this phase to hard 
core dimer coverings of a honeycomb lattice. The critical 
power-law correlations are reduced to exponential corre- 
lations at finite T, but the correlation length may still 
be large. Such a phenomenon would be detectable by 
neutron scattering measurements. 
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APPENDIX A: MEAN-FIELD 
APPROXIMATIONS 

For pedagogical purposes, we examine the zero-field 
TKL Ising model using two mean-field approaches. Com- 
pared with the exact solution (see Fig. [4|), we see that 
such approaches can be quite misleading in the frustrated 
regime J aa /\J a b\ < -I- 
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In the simplest mean-field approximation (MFl), every 
spin is assumed to fluctuate thermally in a mean field de- 
termined by the average magnetizations of its neighbors. 
This leads to a pair of simultaneous equations for the 
spontaneous magnetizations of the a- and 6-sublattices, 
m a = (a a) and mb = (0-5), 



m a = tanh/3(2J aa m a + 2J ab m b ), 
m b = tanh/?(4J ao m Q ). 



(Al) 



where [3 = l/T. Linearizing the tanh function in (er ) 
and (at,) gives the critical temperature, 



T, 



MFl 



J aa 



8 Jab 2 - 



(A2) 



We also present a more advanced mean-field approxi- 
mation (MF2), similar to the method used by Strecka£ 
for the quantum Hciscnbcrg TKL model, in which we 
sum over all eight states of the a-trimers with appropri- 
ate Boltzmann weights instead of treating each a-spin 
independently. The mean-field equations are then 



„4/3J a 



111 



a 



sinh 6/3J ab m b + sinh 2f3J a birib 
3cosh2f3J ab m b ' 



e 4/3J ao cosh 6/3 J a bm 

m b = tanh /3(4 J ah m a ), 



(A3) 



and the critical temperature is given by 



*1F2 



I Jab 



'8(l + 3e 4 -W T e MF2 ) 
( 3 + e 4J a „/rMF 2) • 



(A4) 



T C MF1 and T C MF2 are the dotted and dashed curves, 
respectively, in Fig. [3J The mean-field approximations 
overestimate T c by a factor of two or more. In the regime 
Jaa/\Jab\ ^ 1, T c is dominated by the weak links (J a fe), 
so T c /\ J a b\ tends to a constant, a fact which is captured 
by MF2. However, in the frustrated regime J aa /\Jab\ < 
— 1, both MFl and MF2 predict the wrong behavior of 
T c (see Fig. Also, if an external field is included in 
the analysis, these mean-field approximations predict an 
induced ferromagnetic or fcrrimagnetic moment, but by 
construction, they are unable to capture phases III, IV, 
and the zero-field "In 72" phase in the rich phase diagram 
presented in Sec. IIV1 
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